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The signature for a non-Fermi liquid behavior near a quantum phase tran- 
sition has been observed in thermal and transport properties of many metallic 
systems at low temperatures. In the present work we consider specific exam- 
ples of itinerant ferromagnet as well as antiferromagnet in the limit of vanish- 
ing transition temperature. The temperature variation of spin susceptibility, 
electrical resistivity, specific heat, and NMR relaxation rates at low tem- 
peratures is calculated in the limit of infinite exchange enhancement within 
the frame work of a self consistent spin fluctuation theory. The resulting non- 
Fermi liquid behavior is due to the presence of the low lying critically damped 
spin fluctuations in these systems. The theory presented here gives the lead- 
ing low temperature behavior, as it turns out that the fluctuation correlation 
term is always smaller than the mean fluctuation field term in three as well as 
in two space dimensions. A comparison with illustrative experimental results 
of these properties in some typical systems has been done. Finally we make 
some remarks on the effect of disorder in these systems. 

I. INTRODUCTION 

The description of electronic contribution to the low temperature behavior of metals in 
terms of the Fermi liquids has been highly successful. [Q] The low lying excitations of the 
Fermi liquid manifest themselves in various thermodynamic and transport properties, such as 
the specific heat varying as C v = 7T, a temperature independent (Pauli) spin susceptibility, 
X = 2fi B N(ep), where N(ep) is the density of states at the Fermi energy, a temperature 
dependent electrical resistivity varying as Ap ~ AT 2 , and a linearly temperature dependent 
NMR relaxation rate. T{~ ~ T (Korringa). The values of the these coefficients, such as 7, 
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A, etc., however, are material dependent. For some transition metals these are about one 
order of magnitude larger than in normal metals, and in some compounds containing a large 
concentration of rare-earth or actinide elements such as Ce, Yb or U, these values are about 
a thousand times larger, particularly the value of 7 and the zero temperature susceptibility. 

The normal Fermi liquid behavior as mentioned above is understood within the Landau 
phenomenological theory, where the effect of interaction in a Fermi system is expressed in 
terms of a few parameters which renormalize the physical quantities with respect to their 
free Fermi gas values. For example, the modifications in specific heat, spin susceptibility and 
isothermal compressibility are given by, C v /C® — — = l + (F 1 s /3), x/x° — ( m * / m ) I (1 + F$) > 
and k s /k® = (m*/m)/(l + Fq) respectively. (The superscript denotes the free Fermi gas 
values, other notations are standard). |]] The basic reason for the success of the Landau 
theory is the largeness of quasi-particle life time near the Fermi surface, i.e., r _1 ~| e | 2 ^C e, 
where e = (E — Ep)/Ep. From these relations it is clear that for a certain value of the 
Landau parameters (i.e. Fq, and F\), the corresponding quantities become very large, which 
in turn may indicate a neighborhood of certain phase transition. For example, Fq — > — 1 
implies magnetic instability, or Fq — > 00, a condensation. In the present work we consider 
the Fermi system in the vicinity of such a transition and seek an explanation of the non- 
Fermi liquid behavior of certain substances in this regime. || It seems the Fermi liquid 
theory gives indication of the incoming electronic phase transition as the coupling constant 
is changed, but it does not consider the effect of incipient fluctuations in a self-consistent 
manner. 

There are many examples of electronic phase transitions where the coupling constant 
tunes the transition. These are known as the quantum phase transition. For example 
1 — UN(€f) > gives instability towards ferromagnetism, 1 — U\{Q) > 0, gives antifer- 
romagnetic instability corresponding to a wave vector Q, and n 1 / 3 ^ > 0.26 describes the 
metal insulator transition due to Coulomb correlation as suggested by Mott. These are es- 
sentially zero temperature transitions, however, in general, T c << Tp, where Tp is the Fermi 
temperature of the system. In contrast the classical phase transition occurs at finite tem- 
perature and is described by balance in the energy needed (loss) to create disorder with gain 
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in entropy due to disorder such that the Free Energy, F = U — TS, is reduced. One more 
difference is that the statics and dynamics become correlated in quantum phase transition. 
011 This is principally due to non-commutativity of various terms in the Hamiltonian. For 
example, consider the Hubbard model, 

H = J2e k n k + Uj2n^n ih (1) 

k i 

for correlated electrons. Here the kinetic energy and the U terms do not commute. (Other- 
wise the model will be trivial to solve). Technically, this means that one should introduce 
"time" and the Feynmann time ordering in the functional integral for the partition function. 
The order parameter field becomes "time dependent". The time variable thus acts as an 
extra dimension. This leads to a change in the critical behavior. [[| At a first glance it 
seems that the critical behavior would be the same as that of a D + 1 dimensional classical 
system. However, detailed analysis shows that the critical behavior (or the Upper Critical 
Dimension) depends on the dispersion and damping of the order parameter fluctuations. 
The reason is that the spin susceptibility for a ferromagnet above T c is given by, || 

Uq u>+) « (2) 

For free electron gas 7 = 1/12 and 5 = 1/2. At T c , 1 — UN(ep) —* 0, and therefore to ~ q 3 , 
gives the order parameter dispersion. In case of antiferromagnetism, the staggered spin 
susceptibility is given by, 

x(q+«.^« 1 _ w( g ( ; ) v _^ ) (3) 

In this case, u ~ q 2 at the critical point. A dynamical exponent z is introduced, which 
reflects the change in the static critical behavior. In particular the scaling dimension of the 
quartic interaction is given by e = 4 — (d + z) with z = 3 for ferromagnets and z = 2 for anti- 
ferromagnets. || In field theory z — 1, since uj and q are linearly related, and have the same 
scaling form. At present the application of the renormalization group to quantum critical 
phenomenon in particular the correlation of the static and dynamic behavior is a subject 
of intense activity. We refer the reader to @J^J|,|] for detailed discussion. To summarize, 
the vicinity to the phase transition point and the fermionic nature of correlated electronic 
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system undergoing a phase transition change the nature of the phase transition itself as well 
as the Fermi liquid behavior expected in this system. The reason for this behavior is the 
smallness of the transition temperature, T c compared with the Fermi temperature Tp. This 
aspect gets revealed more clearly as T c — > 0. 

To calculate various physical properties, we take specific examples of ferromagnetic and 
antiferromagnetic transition in itinerant electron system in two as well as three dimensions 
in the limit of vanishing transition temperature near the transition temperature. These two 
examples represent two different types of quantum critical behavior. The basic reason is that 
in ferromagnet the order parameter is a conserved quantity, while in the antiferromagnet 
it is not. This difference is reflected in the dispersion of their respective order parameter 
fluctuation as shown in Eqns. (0) and (0). The microscopic calculation is done within 
the self consistent spin fluctuation theory developed earlier by Ramakrishnan and one of 
10 13] among many others [i"4"-|T6||. For details of the spin fluctuation theory we refer 
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the reader to the monograph by Moriya [17]. A brief review is given in |]18| . We first 
briefly review the spin fluctuation theory and then write expressions for spin susceptibility, 
resistivity, specific heat and the nuclear magnetic relaxation rate. Similar expressions for 
staggered susceptibility and other quantities in antiferromagnets are also written. These 
quantities are then calculated in the limit of large exchange enhancement (i.e. in the limit 
of XpX(T = 0) _1 = a(0) — > 0). Though T c = 0, fluctuation effects are observable well above 
T c . The temperature dependence need not be Fermi-Liquid-like because of the low lying 
fluctuation (Bosonic) degrees of freedom. 

II. SPIN FLUCTUATION THEORY 

The basic motivation for constructing the spin fluctuation theory is the largeness of the 
susceptibility (Stoner) enhancement factor l/a(0). In such case a highly paramagnetic sys- 
tem at low temperature can be considered to be in the vicinity of a magnetic transition. The 
temperature variation of various physical quantities is therefore governed by transverse and 
longitudinal spin fluctuations. Though the order parameter vanishes above the transition, 
the effect of fluctuations is observable well above the transition. There are many equivalent 
formulations available on this idea JT7,TB] . A recent one is due to McMullan, |TI| which is 



based on functional integration using Grassmann variables and collective coordinate trans- 
formation. We briefly summarize our approach and then compile results on some physical 
properties. 

Consider the Landau expansion for the free energy F(M,T), in power of the order pa- 
rameter M, viz. 

F(M, T) = F(0, T) + -A(T)M 2 + -BM 4 - EM (4) 

where H is the field conjugate to M. The temperature dependence of various quantities 
in this theory arises due to A(T) and B. For example, the spin susceptibility for the 
paramagnetic phase is given by 

X '\T) = A(T). (5) 

(A(T) and cx(T) have qualitatively the same temperature dependence and differ only by 
some numerical factors, e.g. A(T) = a(T)/2iV(e_p) for ferromagnets, which we ignore and 
identify A(T) with a(T) now onwards.) Similarly the magnetization in the ordered phase 
is, 

M 2 (T) = (6) 

and the equation of state is given by, 

^L = a{T) + BM 2 . (7) 

The expansion coefficients cx(T) and B have been calculated in various approximation 
schemes. In the Ginzburg Landau theory for classical phase transition, a(T) is taken as 
(T — T c ) and B as independent of temperature. This leads to the Curie Weiss law for the 
susceptibility and the well known mean field critical exponents. In the mean field theory of 
itinerant ferromagnet, 

a MF {T) = (l-UN(e F )y (8) 

and B is again a constant. In this case the temperature dependence of physical quantities 
near T c comes from that of integral over density of states through a Sommerfeld expansion. 
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It is weak, of the order of T 2 /Tp, and therefore it does not give a Curie- Weiss form for the 
spin susceptibility. This issue is tackled in the spin fluctuation theory, where a(T) is given 
by, [0,0] 

a SF (T) = a(0) +u 4 (2D T + 3D L ). (9) 

Here a(0) is the susceptibility enhancement factor at T = 0. This includes the mean field 
part o.mf{T) and the zero temperature part of the fluctuation self energy whose finite tem- 
perature part comprises the second term. Here D T and D L are transverse and longitudinal 
spin fluctuation amplitudes obtained by the internal frequency summation in the diagrams 
shown in Fig. (fj) (a,b,c). The main contribution to the temperature variation of various 
physical quantities is governed by these amplitudes. The factor w 4 in the second term is is a 
dimensionless short range four fluctuation coupling constant obtained after integration over 
fast fermionic degree of freedom. 

The above result has been derived microscopically, within the functional integral scheme 
on a model of interacting electrons. We consider Hubbard model as applied to itinerant fer- 
romagnets and for brevity consider only spin degrees of freedom. Applying the Stratanovich- 
Hubbard functional integral transformation the partition function can be written as 



tr/n%ex P [-Ei^r 

q,m q,m 

~ J du lYl £kn k ,a,u - (^") 1/2 XXCmAm ex P0m M ) + h - c -)]] ( 10 ) 



/n^exp[-/5F(^ m )]. (11) 



where £ 9im is the spin fluctuation field of wave vector q and frequency z m (= 2mm/ (3). Also, 
6k is the kinetic energy of the electrons and U denotes a short range inter atomic repulsion. 
Integrating over the electronic degrees of freedom, we have free energy functional F(£ qtm ) 
for interacting spin fluctuations; that is 

q,m 

Parameters of this model, e.g. fluctuation spectrum, fluctuation coupling vertices are de- 
termined by properties of the underlying fermion system. Since these parameters (e.g. the 
Stoner enhancement factor for ferromagnets or the staggered susceptibility for antiferromag- 
nets) are such that spin fluctuations are low lying excitations, this transformation is specially 
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helpful for an analysis of temperature dependent properties of weak itinerant electron fer- 
romagnets and antiferromagnets. The free energy functional F(£ q . m ) then expanded in 
powers of these fluctuation fields up to a quartic term and a self consistent mean fluctuation 
field approximation (quasi harmonic approximation, or the self consistent renormalization 
scheme of Moriya) can be generated. The mean fluctuation field approximation corresponds 
to the diagram (a), (b) and (c) in Fig. (P, and shown in a compact manner in Fig. @ where 
the double wiggle represents the dressed propagator D(q). The details are given in earlier 
papers [|iT| , |i"2"f . One can also estimate corrections due to higher order fluctuation terms. 
Figs. (HD (d) and (e) and (f) represent typical higher order fluctuation correlated terms. 

III. PHYSICAL PROPERTIES NEAR QUANTUM CRITICAL POINT 

A. Spin Susceptibility 

The self consistent equation for the temperature dependence of a(T) is given by Eq. (|9p 
which is written explicitly as 



a 



(T) = a (0) + A E / dun(aj)lm X (q, oo + ) (12) 



where, A is related to U4, n(uj) = {e UJ,T - 1) , is the Bose distribution function, and 

N(e F ) 

xM = a(T)+6q*- i ^> (13) 
is spin susceptibility for the ferromagnetic case, (u and T are written in units of ep and q 
in units of hp. We have taken h = 1 and = 1). Performing the frequency integral, 

a(T) = q(0) + - £ q{Hv) " f " m) (14) 

n q 2 V 



where, 



An interpolation formula for, 



V = ^WT) + &<?)■ (15) 
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which is valid for small as well as large y is useful in calculating the momentum integral. 
For three dimension, 



A finite ct(0) introduces two regions of temperatures. fLl |. For T < ct(0) one gets the 
standard paramagnon theory results; and for a(0) < T < 1 one gets the classical Curie 
Weiss susceptibility, 

X = Xp/«(T) ~ /4/(aa(0) + T).. (18) 

That is like susceptibility of a collection of classical spins. This feature gets revealed more 
clearly if we put a(0) = in the expression for x(T) and solve the equation self consistently. 
In this case the paramagnon regime (T < a(0)) shrinks to zero and and a classical behavior 
is expected down to T = 0. One is then essentially calculating susceptibility of a ferromagnet 
with T c = 0. Since a(0) is taken to be zero and there is only one region of temperature 
T < 1. In this case, typical y < 1, the limiting form is obtained using the form (fr(y) ~ ^ ( 
valid for y << 1). We then find that, 



/a(T)\V2 , 5 xi/2 

*•-[—) arctan (^^T) 



(19) 



where, gr is a thermal cutoff such that y qT ~ 1. For the form of y given by Eq. fllpp, the 
estimate of the cutoff is ~ Tj/6 or ~ T 1 / 3 . The dominating contribution to o>(T) 
comes from the first term which is given by, T 4//3 . However, since 8 is small, the thermal 
cutoff qx is high w g c (the spin fluctuation energy rises only slowly with q). Thus, a(T) 
rises nearly linearly with T. This is the classical spin fluctuation behavior, first pointed out 



for itinerant ferromagnets by Murata and Doniach. |14| Note that we have assumed T < 1, 
i.e. the system is degenerate. Even so, since the characteristic fluctuation energy a(0) is 
zero, the system behaves classically with regard to spin fluctuations. An estimate of the 
size of the second term is obtained by putting a(T) w T 4 / 3 . We then find it to be of the 
order T 1//3 relative to the first term. Since T 1//3 is not very small, it is essential to do a self 
consistent calculation, particularly in the classical regime which is of interest in the present 
calculation. 
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We have calculated at(T) and other properties in two dimension also. For this we consider 
the same approximate form of the spin susceptibility or the fluctuation propagator as in three 
dimension, the effect of dimensionality is considered only through the phase space in the 
momentum integration. The assumption regarding the form of susceptibility function in 



two dimension is in doubt. It is well known 20 that the Lindhardt function from where 



this functional form has been derived has a different analytic form in two dimension. As far 
as the low momentum behavior is concerned the assumption is closer to the reality if 5 is 
considered to be far smaller than its value in three dimension. For the sake of comparison 
we assume the same value of 5 in three as well as in two dimension. 

Following the same procedure as in three dimensions we find in two dimensions a loga- 
rithmic temperature dependence, 

"< r > = tA§]) < 20 > 

Because of the Bose factor l/(exp(co>/T) — 1), the number of thermal (classical) fluctua- 
tions becomes smaller and smaller as T — ► (i.e. as the T c approaches). This reduces the 
phase space for the fluctuation correlations. In the RNG analysis of Hertz [|J and others || 
this requires the introduction of a suitably scaled 'energy' variable as a degree of freedom 
additional to the three momentum variables. In effect the dimensionality increases, and 
the behavior becomes mean-field. We see this explicitly in our procedure of calculating the 
fluctuation correlation correction perturbatively. The terms involving two or more internal 
thermal spin fluctuations are shown in Fig. (Q). These have been calculated in detail earlier. 
| I"T| It turns out that apart from a numerical factor the two internal thermal spin fluctuation 



term has the same temperature dependence as the mean fluctuation field term. However, 
the three internal thermal spin fluctuation term, « T 2 ln^(l/3a(T)^ in 3D ferromagnet. We 
see that this term is of the order of TlnT relative to the simplest non- vanishing contribution. 
The perturbation expansion therefore converges. 

For a finite T c ferromagnet the mean fluctuation field theory is valid outside the crit- 
ical regime. As the critical regime approaches higher order fluctuation correlations be- 
come comparable to the mean fluctuation term. In the present case, in contrast the 
mean fluctuation field term gives the leading critical behavior. The reason is the follow- 
ing. Suppose for ct(0) = 0, a(T) ~ T A . Then the quantum region T <C ct(T) means 
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T 1 - x < 1 and it occurs only if A < 1. This is not possible, and so one always has the 
other classical (Curie - Weiss) region. Here the fluctuation correlation term is of the form 
T 2 ln(l/o<(T)) ~ T 2 ln(l/T) A < T A . If A ~ 1, the correlation term never becomes more 
important than the mean fluctuation field term. 

In case of Antiferromagnets, the formalism is identical. One replaces the Pauli suscepti- 
bility with the staggered susceptibility for non interacting electron system Xo(Q) 5 f° r brevity 
we retain the same notation for the enhancement factor which is defined in the present case 
as a(0) = Xo(Q)/x(Q)- The expansion of the dynamic staggered susceptibility, Xo(Q + q, 
for small q and small uj around the static staggered susceptibility is also written in the form, 

E 

X (Q + q,^)= — ^ • (21) 

Making similar transformations as for the ferromagnetic case, we get, 

a < T > = a(0) + ^WTi2F) (22) 



where 

(a(T) + 5q 2 



2tt7T 

Thus, for 3D antiferromagnets, 



(23) 



The result turns out to be identical to the ferromagnetic case once we consider only the 
z < 1, where the corresponding momentum cut off turns out to be T 1 / 2 , 



a(T) ( 5 

g c -y — arctan^— g c 



(25) 



and in two dimension there is again a logarithmic behavior, 

( 5 A 
25 n U(T) 
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B. Resistivity 



The electrical resistivity for pure transition and rare earth metals is usually calculated 
within a two band model, where the 'conducting' electrons come from s-band while 
the d-electrons contribute to magnetism. The d-band is assumed to be narrow and the d- 
electrons heavy. The conducting s-electrons scatter from the spin fluctuations corresponding 
to d- electrons. The temperature dependent part of the resistivity due to this mechanism 
for a 3d-ferromagnet is given by, |22|--f24f 

' 1 " ' ' ' ' (27) 



p{T) oc — JqdqJ Imx(q, w )ujn(tu)(l + n{oj))duj 
The frequency integral can be performed by first using the identity, 



J^2 



un(u))(l + n{uj)) 



d_ 
df 



(28) 



leading to, 



p(T)^ T J dqq 3 J duo 



dn(uj) 



quo 



dT q 2 (a(T) + 5q 2 ) 2 + oo 2 



dq q A y(j)'(y) 



(29) 



where, <p'{y) = d< ^^ , and <j)(y) and y are given by Eqns. ( |16D and (|I3D respectively. In the 
limit y < 1 the momentum integral gives, 



P(T) 



T 
26 



a{T)_ x Ja{T) + 5q 2 



(30) 



5 "*V a(T) 

With q c ~ T 1 / 3 we recover the well known result Ap ~ t 5 ^ 3 . However, the self consis- 



tency correction changes the power of temperature. Similarly for two dimensions, 

P(T) 



T 
J 



Qc 



la(T) 



arctan 



s 



(31) 
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The case of 3D antiferromagnets formalism is similar except for the power of q in the 
momentum integral. This is due to the fact that the small momentum expansion is not done 
around q = but around q = Q, the antiferromagnetic wave vector. The result is, [0 



p(T) oc — J q 2 dq J imy(Q + q, uj + )um{uj){l + n(uj))duj. 



(32) 
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Following the same steps as for the ferromagnetic case, we get, 

;i + i2z) 



P (T)<xr/rtii 



2z(l + 6z) 2 

where z is given by Eq. (|23|) The result in the limit of z < 1 is, 



(33) 



P(T) = | 



/a(T) 
q c — \l — - — arctan 



5 



S 



a(T) 



(34) 



Similarly for two dimensions, p(T) = ^ln 



6q? 



26 \ a(T) I ■ 



C. Specific heat 



The spin fluctuation contribution to the free energy within the mean fluctuation field 
approximation (or quasi harmonic approximation) is given by, fl3| 



3T 

Afi = -y E Ml - + AT £ D q , m ,}. 



(35) 



Where -D 9im is the fluctuation propagator which is related to inverse dynamical susceptibility, 
and Xqm i s the f ree Fermi gas (Lindhardt) response function. The argument of the logarithm 
is related to inverse dynamic susceptibility Considering only the thermal part of the integral 
and ignoring the zero point part, we perform the frequency summation and obtain, 
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^^Thermal ~ — / . I , , -, 

7T ^ Jo e"/ 1 - 1 



dw , nuj/Aq 

arctan W) + ^ } ' 



(36) 



Integrating over frequency, we get, 



^Thermal = 3T E ( ln ~ (V ~ ^ +V ~^ 



(37) 



where, y is given by Eq. (|15|) . Once the free energy correction is known, the specific heat 
correction is given by 



AC, 

kn 



-T 



d 2 Atl 



dT 2 

-3T 2 £ 



,2 dy 9 2 y dy 2 d<p(y) 

^TdT dT 2)<P[y) { 8T ) dy 



6 / q 2 dq{<P'(y)( 



7r 2 7 dT 



g d 2 a{T) 
7r 2 7 dT 2 



}• 



(38) 
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Making the small y approximation and introducing the appropriate cutoff, 



AC, 



2ir 2 



dT 2 

\df~) Jo 



dq 



dq- 



a{T) + 5q 2 
qc 



a(T) + 5q 2 ) 2 



o 



dq q l 



,d 2 a(T) t nrT1 da(T)\ 
~df~) 



- 1AT2 -^ + 2T: 



la(T) 



arctan I 



5 



+ T 2 (da(T)\ 2 



25 V dT 



1 



a{T)5 



arctan 



a(T) 



Q< 



W a(T) 



a(T) + 5q 2 



3 ' 



(39) 



The last result is obtained after the momentum integration. Approximately the terms can 
be arranged as, 



2tt 2 



ql + T 



V dT J 



1 



T^T) 



T 



dT 2 



+ 2T- 



dT ) 



Qc 



(40) 



The first term gives the classical result (for constant cutoff) , the second dominant term gives 
leading temperature correction, the last term is about two order of magnitude small at the 
temperature range of interest. Similarly in two dimensions, 



25 V dT 2 

T 2 ( da(T) \- 
~25a{T)\ dT ) 



1 + 



dT 
a(T) 



In 



a(T) 



+ 



(41) 



The calculation for antiferromagnet is identical except that y is replaced by z in Eq. ( j40|) 
and the final equation in terms of a(T) turns out to be identical except that the temperature 
dependence of a(T) is different in an antiferromagnet. 



D. Nuclear Spin Relaxation Rate 



The nuclear spin relaxation rate in metals is given by the Korringa relation, p 



^ « (^y (42) 

T{T \ H J K ' 

which essentially tells that l/Ti is proportional to the square of the static spin susceptibility 
of metals, which in turn, is independent of temperature for most normal metals. However, 
it has been pointed out by Moriya p6[ long ago that this relation gets modified in presence 



of electron correlations. The nuclear spin-lattice relaxation rate in metals is given by, 
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1 „ Imx + (q,uJ 



TlT E^f^ (43) 

where ujq is the nuclear magnetic resonance frequency which is taken to be very small (— > 0) 
in the problem of nuclear spin relaxation rate. Substituting the expression for X~ + { ( l7 UJ o) 
and taking the limit, we have, 

T\T ~ ^ q(a(T) + 5q 2 ) 2 ' (44) 

for a ferromagnet in three dimension. After the momentum integration, the result is 
(TiT) -1 ~ a(T) _1 . For normal Fermi liquid a(T) is constant but in the present case it 
varies as T 4 ' 3 . This leads to a non- Fermi liquid behavior again. Similar calculation is done 
for antiferromagnets. 

E. Effect of Disorder 

The effect disorder can be included in the above mentioned formalism by modifying the 
propagators and vertices in diagrams for the spin fluctuation self energy. This has been 



done in our earlier papers. [27| , 28) In presence of disorder the electron moves in random 
way getting scattered from impurities repeatedly. This introduces a finite mean free path 
for electron, and a finite life time r in the electron propagator, which also modifies the 
free particle-hole propagator (diffuson), free particle-particle propagator (Cooperon) and 
electron-spin fluctuation vertex. The correction to ot(T) to leading order is given by, 

a{T)=a SF {T)-a d (T), (45) 

where ad(T) is correction due to diffusive modes. It is given to the leading order in l/e F r 
as 

a d (T) ~ [1 - sJ(27rrT)}/(e F r) 2 , (46) 

for T C 1/r and it vanishes otherwise. Clearly the disorder introduces a new energy scale 
(1/t), in the lowest temperature range as shown schematically in Fig. (|3|). In the case of 
non- vanishing a(0), in the case of ferromagnet in three dimension, the susceptibility inverse 
a{T) behaves as ~ T 1 ' 2 for T < 1/r, as T 2 /a(0) for 1/r < T < a(0), and as T 4 / 3 for 
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ct(0) < T < 1. Similarly the resistivity correction Ap(T) behaves as ~ T 1 / 2 for T < 1/r, as 
T 2 / ya(0) for 1/r < T < a(0), and as T 5//3 for a(0) < T < 1. In the case of a zero-T c system 
the quantum fluctuation regime (T < a(0)) vanishes and the other two regimes merge. At 
lowest temperature the effect of diffusive mode seems to give the dominant contribution, i.e. 
a(T) ~ T 1 / 2 and Ap(T) ~ T 1 / 2 but a more detailed analysis is needed. 

F. Summary 

The following table summarizes our results. In the first column Fermi liquid theory 
results are written, other columns compile the fluctuation theory results. These results 
are presented in three rows for each property. The first row gives results from a non-self- 
consistent calculation, i.e., for example, when only the first term in Eq. (|19|) for a(T) is 
considered but with a proper momentum cut off. This behavior is expected in the extreme 
low temperature range. These results are in general known, but presented in the coherent 
form for the first time. The second row gives these results with the temperature dependence 
of oi(T) taken in account and the integration performed with the functional form for <f>(y) 
valid for all y but approximated by Eq. fllBp. The power of temperatures so obtained 
depends slightly on the temperature regime considered (i.e. whether T is in 10~ 3 — 10 -2 
or otherwise). The third row gives the classical spin fluctuation results, where the Bose 
factor n(uj) is approximated as T/u, (effectively the first row with a constant cutoff). The 
experimental results are expected to lie between those given in rows one and three. 
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IV. EXPERIMENTAL RESULTS: 



In this section we give examples of materials exhibiting non-Fermi liquid behavior at low 
temperatures and also compare some results with a theory presented above. The most popu- 
lar example of system showing non-Fermi liquid behavior is, of course, the high temperature 
superconductors [pj . It seems, however, the effective low dimensionality, the specific nature 
of the density of states and structural aspect of the Fermi surface (nesting etc.) play im- 
portant role in this system. We therefore want to consider examples from three dimensional 
correlated electronic system in the neighborhood of electronic phase transition. 

The next example is that of phosphorus-doped silicon (Si:P). [30| This material goes 
through a Mott insulator to metal transition as the doping by P increases. At n c ~ 3.7 10 18 
P per cm 3 there is a metallic state. The spin susceptibility of Si:P gets enhanced and 
becomes strongly temperature dependent as the metal-insulator transition is approached 
from the metallic side. The T-dependence observed does not fit the 1/T behavior expected 
for weakly interacting localized spins either. Moreover, the spin-lattice relaxation times in 
barely metallic Si:P is strongly temperature dependent, T^ 1 ~ T -1 / 2 , similarly the correction 
to the zero-T conductivity, a ~ T 1 / 2 , in this material. There are theories which associate 
these anomalies to spin fluctuations induced due to incipient localization. There is a subtle 
interplay of disorder and correlation effect in this material. Only spin fluctuation kind of 
theory will not work. 

Some transition metal (and also some actinide) inter-metallic compounds show a low 
saturation moment per transition metal atom and a low magnetic transition temperature T c 
compared with conventional ferromagnetic materials like Fe, Co and Ni. These compounds 
are known as weak itinerant electron ferromagnets. The prototype examples are, ZrZn 2 , 



Ni 3 Al, and Sc3in. [18| Their low temperature properties have been discussed within spin 



fluctuation theories for a long time. [|17] Here we compare the specific heat behavior of Sc3in 
above T c with our present calculation (Fig. (f|)). The experimental curves are due to Ikeda 
and show a good fit to the theory with Ap(T) ~ T, and C V {T)/T ~ (T - T c )-°' 25 . 
The example of MnSi is interesting from the perspective of the present work. The 
material has a transition temperature around 30A". As the hydrostatic pressure is applied 
the T c decreases continuously and collapses towards absolute zero at p c = 14.6kbar. This 
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is an example where an approach to a zero temperature quantum phase transition can be 
observed as a function of pressure. This has been done by the Cambridge group. The 
deviation of the resistivity curve from the T 2 behavior gets pronounced as p c is approached. 
In Fig. (H)we have compared Ap/T 2 with a power law temperature dependence as suggested 
in our calculation. 

We have compared the nuclear spin relaxation rate of 27 Al in Ni^Al as a function of 
temperature. |£J This material has a transition temperature about 41 K and shows all 
other characteristic properties of weak itinerant ferromagnet. |33] The low field data fit to 
power law T -0,89 for {T{T)~ l as shown in Fig. 

The heavy fermion material CeCug is non-magnetic. On alloying with Au the lattice 
expands and an antiferromagnetic order is observed in CeC\x§- x h\i x above a critical con- 
centration x c ~ 0.1. The Neel temperature of the anti-ferromagnetic heavy-fermion alloy 
CeCu5.7Auo.3 can be continuously tuned to zero with increasing hydrostatic pressure. At the 
critical pressure the specific heat has been fitted to C/T ~ In T /T curve. [P5[] We analyze 
the data again and fit the curve to our prediction (T 58 corresponding to the temperature 
range of interest) in the Fig. fl7|). 



V. CONCLUSION 



We have calculated the temperature dependence of various physical properties near the 
quantum phase transition point. The results hold for electronic phase transitions with a 
finite T c also. This is clear, as the results from the first rows in Table (1) match with some 
well known results in literature. However, they are pronounced and a clear non-Fermi liquid 
behavior is obtained when T c — > 0. Our results are perturbative, but as discussed in the text 
the fluctuation correlation term is always smaller than the mean fluctuation field term. The 
behavior of these quantities is different in ferromagnet from the antiferromagnetic system. 
This is a reflection of the fact that the order parameter fluctuations have different form 
of dispersion in these systems. Finally we made some remarks about the inclusion of the 
effect of disorder near quantum critical point within the spin fluctuation formalism. The 
present approach can be applied to other systems also. One only needs an appropriate form 
of the order parameter correlation function to calculate various quantities. For example this 
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approach can be applied to systems with a pseudo gap in the excitation spectrum and 
also with phonon like dispersion as it happens in short coherence length superconductors. 
For example in 2D short coherence length superconductors it has been shown through Monte 
Carlo simulations that the relaxation rate varies as the spin susceptibility, |37| (a(T)) -1 of 
the system, which matches with our result on 2D antiferromagnets. 
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FIGURE CAPTIONS 



FIG. 1. Self energy diagrams for the spin fluctuation propagator. 
FIG. 2. Self energy in the mean fluctuation field approximation 

FIG. 3. Schematic diagram of the various temperature scales involved in the disordered material 

FIG. 4. Plot of C v /T as a function of T — T c for Sc^In. The experimental points are represented 
by circles and the solid line represents the theoretical fit. 5 = 1/12, 7 = 1/2, Tp ~ 1000K. 

FIG. 5. Plot of p(T) as a function of T for MnSi. The experimental points are represented by 
circles and the solid line represents the theoretical fit. S = 1/12, 7 = 1/2, Tp « lOOO-fT. 

FIG. 6. Plot of (T\T)~ l as a function of T for Ni%Al. The experimental points are represented 
by circles and the solid line represents the theoretical fit. 5 = 1/12, 7 = 1/2, Tp w 1000K . 

FIG. 7. Plot of C v as a function of T for CeCu^,jAuQ,3 at a pressure of 8.2 kbar. The 
experimental points are represented by circles and the solid line represents the theoretical fit. 
5 = l/2vr,7 = l,Tp w 5K. 

TABLE CAPTIONS 

TABLE I. Summary of the temperature dependence of various thermal and transport properties 
near a quantum phase transition point. The first row gives the non-selfconsistent calculation 
scheme. For ferromagnets the upper cutoff for q is T 1 / 3 while for antiferromagnets it is T 1 / 2 . 
The second row gives the selfconsistent calculation scheme results. The upper cutoff for q has 
been taken to be 1. The range of temperatures in which these exponents have been calculated is 
T = 10~ 3 to 1CT 2 . The third row gives the classical spin fluctuation results, (i.e. the first row 
with qx as a constant.) 
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